Probing the dark matter profile of hot clusters 
and the M — T relation with XNLNL- Newton 



E. Pointecouteau, a M. Arnaud, a G.W. Pratt b 

l SAp/CEA, I'Orme des Merisiers, 91191 Gif-sur-Yvette, France 
h MPE, Giessenbachstrafie, 85748 Garching, Germany 



Abstract 

We present results based on XMM-Newton observations of a small sample of hot 
galaxy clusters. Making a full use of XMM.-Newtori's spectro-imaging capabilities, 
we have extracted the radial temperature profile and gas density profile, and with 
this information, calculated the total mass profile of each cluster (under the assump- 
tion of hydrostatic equilibrium and spherical symmetry) . Comparing the individual 
scaled total mass profiles, we have probed the Universality of rich cluster mass pro- 
files over a wide range of radii (from 0.01 to 0.7 the virial radius). We have also 
tested the shape of cluster mass profiles by comparing with the predicted profiles 
from numerical simulations of hierarchical structure formation. We also derived the 
local mass-temperature (M — T) scaling relation over a range of temperature going 
from 4 to 9 keV, that we compare with theoretical predictions. 

Key words: Clusters of Galaxies, Large Scale Structure, X-rays, 
XMM-TVeioton Observations 



1 Introduction 



In a framework where cosmic structure forms under the effect of gravitation 
alone, the population of galaxy clusters should be strongly similar. A simi- 
larity of shape is expected for the dark matter distribution in clusters due to 
the process of gravitational collapse. The gas "follows" the distribution of the 
dominant dark matter component within the gravitational well, and is in ap- 
proximate hydrostatic equilibrium (between major mergers). Then the matter 
in clusters (be it baryonic or dark) has a similar distribution from cluster to 
cluster. Scaling laws are another aspect of cluster similarity and are a direct 
prediction from scenarios of structure formation based solely on gravitation. 
They express the correlation between physical quantities (such as temperature, 
luminosity, entropy,...), and in fact make clusters a two parameter population 
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which depends only on the redshift and the total mass of each object. Each 
single other quantity can then be expressed as follows: Q(z) = a(z)M a . 



The properties of cluster similarity are a powerful statistical tool for the study 
of the cluster population, and a ri ch source of informat ion on structure forma- 
tion and evolution. (For details see lArnaud et al.l (J2004J) and reference therein). 



The paper presents results in the currently-favored concordance cosmological 
model: H = 70 km/s/Mpc, Vt m = 0.3 and fi A = 0.7. 



2 Universality of the mass distribution 



We have investigated the shape of the dark matter profile in clusters of galaxies 
using a sample of five rich, relaxed nearby clusters with temperatures ranging 
between 4 and 9 keV. The different steps in the computation of each mass 
profile are: 

• The computation of an azimuthally-averaged surface brightness (SB) profile, 
which is then parameterize d with either (a) the sum o f n /3-models (n = 



2 or 3, BB or BBB model, iPointecouteau et al.l 120041 ) or (b) a modified 



double /3-model allowin g a peaked density dis tribution toward the cluster 
center (i.e KBB model. IPratt fc Arnaud1l2002h . 

Spatially-resolved spectroscopic fits of annular regions to obtain a temper- 
ature profile. The output projected temperature profile is then corrected 
for projection effects and for the effects of PSF blurring. This is achieved 
through a Monte Carlo simulation in which the projected profile is ran- 
domly sampled 1000 ti mes within the obs erved errors, parameterized using 
the functional form of lAUen et al.l (|200ll ). and subsequently corrected for 
projection and PSF effects. 

Under the assumptions of spherical symmetry and hydrostatic equilibrium, 
these density and the temperature profiles can then be used to derive the 
total mass p rofile. More deta il ed explanat i ons o f thos e different steps can 
be found in IPratt fc Arnlul (|2002l I2003L l2005h and IPointecouteau et al.1 
f|2004h . 



We have used our derived mass profiles to test different models for the dark 
matter distribution . We tested an isothe rmal sphere model, a MQGSL profile 
(jMoore et allll99flh and a NFW profile (|Navarro et al.l . ll997t ). The NFW pro- 
file proved to be the best representation of our current data. We adopted this 
as the parametric description for our cluster mass profiles. In fitting the mass 
profiles, the free parameters were the mass at a given overdensity with respect 
to the critical density of the Universe, M(S), and either the scale radius, R s , 
or the concentration parameter, c. This approach allows us to work directly 
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Fig. 1. Scaled mass profiles. Each profile is scaled with respect to its virial radius 
-R200 an d virial mass M200 ( see text). The solid black line corresponds to the average 
NFW profile and the two dashed lines are the associated standard deviation. 
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Fig. 2. The M — T relation at 5 = 2500 from 5 hot relaxed clusters observed with 
XMM-iVewfon (red filled dots with associated error bars); the associated best fit is 
overplotted as a solid red line. The dashed black line is the pre diction from numerica l 
simulations based purely on gravitationaljstructure formation lEvrard. et al. (1996). 



Table 1 

Details on the sample 



Name 


z 


Rev. 


^exp 


n e (r) model 


A478 


0.088 


401,411 


48/41/37 


BBB 


PKS0745 


0.103 


164 


10/10/- 


KBB 


A1068 


0.137 


633 


19/20/15 


BBB 


A1413 


0.143 


182 


24/25/10 


KBB 


A2204 


0.152 


322 


20/20/13 


KBB 



Columns: (1) Name. (2) redshift. (3) XMM-Newton revolution. (4) Effective ex- 
posure time for EMOS1/EMOS2/EPN. (4) Parametric model used for the density 
profile modelisation. 

with the mass at a chosen overdensity. Each pair of parameters led to identical 
results. 

We then scaled each total mass profile according to the cluster virial radius 
(defined as the radius where the density is equal to 200 times the critical 
density of the Universe at the cluster redshift, -R200) an d the virial mass 
(the mass at this radius, M200). Our small sample includes the following 
clusters: A478, PKS0745, A1068 , A1413 and A2204. The result s on A 1413 
and A478 have been respe ctively published by Pratt fc Arnaudl (2002) and 



Poi ntecouteau et al.l ([20041 ) . PKS0745 was originally analyzed bv lChen et al. . 

but we reprocessed the whole data set checking that our result are 
consistent with the previously published ones. Nothing on A1068 and A2204 
has yet been published. The details on clusters are gathered in Table 1). The 
five scaled mass profiles are plotted together in Fig. 1. The five profiles scale 
very well over a large range of radii from about 0.01 to 0.7 the virial radius, 
exhibiting a strongly similar shape as expected for an underlying universal 
matter profile in clusters. 



3 The M - T relation as seen by XMM- Newton 



In order to further check the properties of cluster similarity with our sam- 
ple, we have investigated the mass to temperature (M — T) relation. For 
easy comparison with previous work, we focus on an overdensity of 5 = 2500 
with respect to the critical density of the Universe. The "virial" temperatures 
used to determine the M — T relation were obtained from direct spectro- 
scopic fits over a region from r > 0.1R V , chosen to avoid the cooling cores 
found in many of these clusters. We then performed a linear regression fit 
in the log M — log kT plane, taking into account both the errors in M and 
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kT. With only five XMM- Newton clusters we are able to derive the follow- 
ing M — T relation: M 2500 = 1.45 x 10 13 (A;T/keV) L56±a16 M (see Fig. 2). 
This result is compatible in terms of slope with the slope of 1.5 predicted 
by numeri cal simulations relyin g on gravitation-based scenarios of structure 
formation ( Evrard. et all Il996l ). However, the normalization of the observed 
relation is lower by a factor of about 30%, a problem encountered in several 
previous works jNevalainen et al. . 2000l : Finoguenov et al. . 2001 : Allen et al. 



2001 : Sanderson et all 120031 ). We then added the published Chandra dat 



a on 



four rich (kT > 5 keV) clust ers not already in our sample (A2390, A1835, 
RXJ1347-1145, MS2137-2353; lAUen et all bnOlh and re-derived the M — T 
relation, finding M 250 o = 1-46 x 10 13 (fcT/keV) L55±a15 M . The fit is almost 
unchanged, perhaps due to the larger uncertainties of the Chandra measure- 
ments with respect to the very precise XMM-Newton measurements. 



Thus at 5 = 2500, our derived M — T relation is compatible in terms of slope 
with the prediction from scenarios of structure formation based on gravitation 
alone. 



4 Conclusion 



With this work we have used XMM- Newton's excellent imaging and spectro- 
scopic capabilities to observe clusters of galaxies up to about half of the virial 
radius (r ~ 0.5i? 20 o)- In this study we were able to: 

(1) probe the universality of the dark matter distribution in a sample of 5 
rich, nearby galaxy clusters from 0.01 to 0.7 times the virial radius. The 
five scaled mass profiles present a very similar shape which is fairly well 
described by an NFW density distribution. 

(2) investigate the M—T relation at the overdensity of 2500 times the critical 
density of the Universe over a range of temperature going from about 
4 keV to 9 keV. The derived M — T relation has a slope compatible 
with the expected slope of 1.5 from the numerical simulations based on 
gravitational heating. The precision of the XMM-Newton measurements 
allows us to put tight constraints on the M — T slope with only 5 clusters. 
Still our observed normalization for this relation is lower by a factor of 
about 30% compared to theoretical predictions. 

Since the mass is not measured directly, but is derived from the density and 
the temperature profiles, a key issue for this kind of work lies in the accuracy of 
the determination of the true (meaning deprojected and PSF corrected) tem- 
perature and surface brightness profiles. The shape of Chandra temperature 
profiles (Vikhlinin, this meeting) and the XMM-Newton profiles we derived 
from our observations have therefore to be cross-compared. This is essential 
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to assess the reliability of temperature profiles in mass profile determination 
and in the calibration of scaling laws. In a forthcoming paper, we extend our 
sample coverage to lower temperature (Pointecouteau, Arnaud & Pratt, in 
preparation). There we will investigate in greater detail the properties of the 
scaled profiles, the slope and normalization of the M — T relation, and com- 
pare our results to recent numerical simulations involving non-gravitational 
processes. 

The quality of current X-ray data makes it possible to calibrate the scaling 
laws in the local universe (the highest redshift cluster of our sample is A2204 
with z = 0.1523) with a high precision which is fundamental for: (i) further 
studies concerning the evolution of scaling laws, and therefore for the studies 
of structure formation and evolution; (ii) surveys at different wavelengths like 
blind SZ surveys, which will make use of calibrated X-ray scaling laws as a 
bridge to calibrate SZ scaling laws such as the Y — M relation (roughly the 
SZ flux to total mass relation), needed for cosmo logical studies based on SZ 
statistics. 

Indeed for such kind of studies, biases and/or systematic errors could be in- 
duced by the thermodynamical state of the clusters, or by some hypothesis 
or processes in the data analysis. Such kind of biases and systematics have 
to be carefully investigated and quantified to assess the calibration of scaling 
laws, and to make their use reliable on large cluster samples for cosmological 
studies and structure formation studies. 
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